arXiv:astro-ph/0512355v2 25 Jan 2006 


Exponential Cardassian Universe 


Dao-jun Liu, Chang-bo Sun, Xin-zhou Li 

Shanghai United Center for Astrophysics(SUCA), 

Shanghai Normal University, 100 Guilin Road, Shanghai 200234, China 


Abstract 

The expectation of explaining cosmological observations without requiring new en¬ 
ergy sources is forsooth worthy of investigation. In this letter, a new kind of Car¬ 
dassian models, called exponential Cardassian models, for the late-time universe are 
investigated in the context of the spatially flat FRW universe scenario. We fit the 
exponential Cardassian models to current type la supernovae data and hnd they 
are consistent with the observations. Furthermore, we point out that the equation- 
of-state parameter for the effective dark fluid component in exponential Cardassian 
models can naturally cross the cosmological constant divide w = —1 that observa¬ 
tions favor mildly without introducing exotic material that destroy the weak energy 
condition. 
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1 Introduction 


The current accelerating expansion of the universe indicated by the astro¬ 
nomical measurements from Supernovae type la (SNela) [1] (see [2,3] for most 
recent results) as well as accordance with other observations such as the cos¬ 
mic microwave background (CMB) [4] and galaxy power spectra [5] becomes 
one of the biggest puzzles in the research held of cosmology. There are lots of 
approaches to unriddle this puzzle. One popular theoretical explanation ap¬ 
proach is to assume that there exists a mysterious energy component, dubbed 
dark energy, with negative pressure, or equation of state with w = p/p < 0 that 
currently dominates the dynamics of the universe (for a review see [6]). Such 
a component makes up 70% of the energy density of the universe yet remains 
elusive in the context of general relativity and the standard model of particle 

Email address: kychz@shnu.edu.cn (Xin-zhou Li). 


Preprint submitted to Elsevier Science 


5 February 2008 





physics. In recent years, many candidates for dark energy have been explored. 
Besides a cosmological constant, one popular candidate source of this missing 
energy component is a slowly evolving and spatially homogeneous scalar held, 
referred to as ”quintessence” with tc > —1 [7] and’’phantom” withtc < —1 [8], 
respectively. Since current observational constraint on the equation of state of 
dark energy lies in a relatively large range around the so-called cosmological 
constant divide wx = —1, it is still too early to rule out any of the above 
candidates. 


On the other hand, general relativity (GR) is very well examined in the solar 
system, in observation of the period of the binary pulsar, and in the early 
universe, via primordial nucleosynthesis. However, no one has so far tested in 
the ultra-large length scales and low curvatures characteristic of the Hubble 
radius today. Therefore, it is a priori believable that Friedmann equation is 
modihed in the very far infrared, in such a way that the universe begins to 
accelerate at late time. Freese and Lewis [9] construct so-called Cardassian 
universe models that incarnates this hope. In Cardassian models [9,10,11], the 
universe is flat and accelerating, and yet contains only matter (baryonic or not) 
and radiation. But the usual Friedmann equation governing the expansion of 
the universe is modihed to be 


= 




( 1 ) 


where p consists only of matter and radiation, H is the Hubble ’’parameter” 
which is a function of time, a is the scale factor of the universe, and G = l/rripi 
is the Newtonian gravitational constant. Note that as required by inhation sce¬ 
nario and observations of CMB, the geometry of the universe is hat, therefore, 
there are no curvature terms in the above equation. Perhaps the most inter¬ 
esting feature of Cardassian models is that although being matter dominated, 
they may be accelerating and may still reconcile the indications for a hat uni¬ 
verse (Ct = 1) from CMB observations with clustering estimates that point 
consistently to Cm,o = 0.3 with no need to invoke either a new dark component 
or a curvature term. The expectation of explaining cosmological observations 
without requiring new energy sources or certainly worthy of investigation. 

For any suitable Cardassian model, there exist at least three requirements 
that should be satished. Firstly, the function g{p) should returns to the usual 
form p at early epochs in order to recover the thermal history of the standard 
cosmological model and the scenario for the formation of large scale structure. 
Secondly, g{p) should takes a diherent form at late times 2 : ~ 0{1) in order 
to drive an accelerated expansion as indicated by the observation of SNela 
[1,2,3]. Finally, the classical solution of the expansion should be stable, i.e. 
, the sound speed of classical perturbations of the total cosmological fluid 
around homogeneous FRW solutions cannot be negative. 
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For the original power-law Cardassian model g{p) = p+Bp^, where B and n < 
2/3 are two constants [9], the second term in the above equation behave like 
an effective cosmic dark energy component that drive the universe accelerate 
at late times and is negligible at early epochs. However, the sound speed of 
cosmological fluid in this model is not guaranteed to be positive. So this model 
should only be considered as an effective description at scales where the sound 
speed is positive [14]. The generalized power-law Cardassian models, such as 
MP Cardassian model [11], satishes all of above requirements, resorting to an 
additional parameter. 

In this paper, we investigate another kind of Cardassian models, dubbed Ex¬ 
ponential Cardassian models. In the next section, we construct two models 
that embodies the elegant idea of Cardassian universe. In section 3, we £t the 
parameters of the models to the type la supernovae observations. The issue 
on the sound speed of the cosmological fluid is addressed in section 4. And 
in the last section, a brief discussion is included. 


2 Exponential Cardassian models 


2.1 Model I: a simple exponential Cardassian model 


In this subsection we consider the following simple version of the Exponential 
Cardassian model; 





SttG 

pexp 



( 2 ) 


where pcard is a characteristic constant energy density and n is a dimensionless 
constant. For the late-time evolution of the universe we neglect the contribu¬ 
tion of radiation. The current acceleration expansion of the universe requires 
that 


3n 


) n 

’ 


(3) 


where po is the current value of energy density of matter p in the universe, 
which keeps conserved during the expansion of the universe, i.e. 

p + 3H{p + p) = 0. (4) 
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Therefore, the evolution of matter takes the ordinary manner 


p = Po{l + zf. 


( 5 ) 


Obviously, at early times, p is much larger than the characteristic energy 
density pcardi g{p) P, i-e. Eq.(2) recovers the standard Friedmann equation. 
In terms of Eqs.(2) and (4), the effective pressure of total fluid pr takes the 
following form: 


Pt = P 


dgjp) 

dp 


g(p)- 


( 6 ) 


In the light of the observation that depend only on the scale factor, in the 
late time regime of the universe hlled with just matter, Cardassian models 
undifferentiated from the effective dark fluid models with the equation of state 
Px = wxpx, where Px = g{p) - P and px = Pt- 

Using the dynamics of p, we change the Eq.(2) into the form 


= Hi 


Om,o(l + z)'^ + (1 ~ ^m,o)fx{z) 


(7) 


where the current value of Hubble parameter Hq is often denoted as 100 h km s ^Mpc ^ 
in which h is a dimensionless constant and 


r / \ T z') 

Jx[z) - 


1 -n 


m,0 






d- 


( 8 ) 


It is easy to check that fx{z = 0) = 1. Therefore, in this model, besides the 
density parameter Qrn,o, there is only one free parameter n for the effective 
dark fluid as in the original power-law Cardassian model [9]. 

Using Eq.(6), we can get the effective parameter of equation of state wt for 
the total cosmological fluid in this model 


WT = ^44- = = nln(fI^,o)(l + 

g{p) p 


(9) 


However, the equation-of-state parameter of effective dark fluid wx becomes 


Pxjp) _ Pt{p) _ nln(fI^^o)(l + 2:) 

Px(p) g{p)-p l-exp[lnH,„,o(l + 2:)"^’"!’ 


( 10 ) 
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Fig. 1. The evolution of if as a function of redshift z in the context of exponential 
Cardassian model where we have choose the parameters n = 0.75, Qm,o = 0.40. 
It is obvious that the equation-of-state parameter of effective dark fluid wx (down 
panel) cross the cosmological constant divide w\ = — 1, whereas that of total ffuid 
is greater than —1 (up panel) up to now. 

In Fig.l, we show the evolutions of the equation-of-state parameter for the 
total cosmological fluid wt and that for the effective dark fluid. Clearly, 
Wx (down panel) cross the cosmological constant divide wa = — 1 recently, 
whereas wt is always greater than —1 (up panel) up to now. However, we 
note that both wt and wx will be less than < — 1 and trend to negative inhn- 
ity with the expansion of the universe. Therefore, the universe will suffer from 
the fate of so-call ” big-rip” in the future [15]. Interestingly, wx is a negative 
constant in the original power-law Cardassian model [9]. Therefore, this can 
lead to discrepancies between exponential and power-law models. 


2.2 Model II: a modified version of exponential Cardassian model 


We now consider the following modihed version of the original Cardassian 
model: 


dip) 


(P + Pcard) exp 


QPcard 
P pcard 


n-\ 




( 11 ) 


where Pcard is a characteristic constant energy density and q and n are two 
dimensionless positive constants. Obviously, at early times, p is much larger 
than the characteristic energy density Pcard) g{p) P, i-e. Eq.(ll) recovers 
the standard Friedmann equation. In terms of Eqs.(ll) and (6), the current 
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acceleration expansion of the universe requires that 


o-Q 


1 -3n 


1 + o-Q 



( 12 ) 


where we dehne the dimensionless quantity a = pjPcard = Po/Pcardi)- + 2 :)^ = 
(To(l + 2 :)^.According to the assumption that the universe is flat, there is a 
relationship between the four parameters o'o, 0. and n that 


(1 + (To) exp 


1 + CTo 


n-\ 


(Toff 


-1 

m,0* 


(13) 


Using the dynamics of p and Eq.(13), the function fx{Z) in Eq.(7) that de¬ 
scribe the evolution of the effective dark fluid becomes 


fx{z) = 


Q 


m,0 


n 


m,0 


<^0 ^ + (1 + 
-(1 + ^)'}. 


exp 


-(1 -|- (To)"' ln[(l -|- (To )ffm,0 


[l + ao(l + ;^)3]' 


(14) 


Combining Eqs.(14) and (13), it is easy to check that fx{z = 0) = 1. There¬ 
fore, in this modihed model, if we assume a priori value for ffm.o 1 there is 
only two free parameters q and n for the effective dark fluid as in the MP 
Cardassian model [11]. 

Using Eq.(6), we can get the effective parameter of equation of state wt for 
the total cosmological fluid in this model 


Wt = —1 -|- 


(T 


1 + a 


1 + n 


1 + CTp 

I + a 


ln[(l -|- (Jq 


(15) 


and the equation-of-state parameter of effective dark fluid wx becomes 


—l + an\ 


1 ln[(l-7ao^)fim,o. 

1 + a 

1 — exp 


rin[(l -7ao^)flm,o] 

)] 


(16) 


In the limit of (T ^ 0, both wt and wx trend to —1, therefore, the expansion 
of the universe will speedup forever and in the end the universe becomes de- 
Sitter universe asymptotically. For the current epoch a{z = 0) = (To, therefore, 
the current value of wx is determined by 


wx{z = 0) 


-1 -I- nap ln[(l (Tq ^)nn^,o[ 
(1 -|- (To)(l — flm,o) 


(17) 
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Fig. 2. The equation-of-state parameters wt and wx as two functions of a in the 
modified version of exponential Cardassian model for different values of parameter 
q, where we have chosen ^2^,0 = 0.30, n = 5. 

Obviously, the value of wx{z = 0) can be greater or less than —1. Furthermore, 
from Fig.2, we can see the evolution of the wx{z). It is interesting that, for 
some values of parameters, wx can cross the cosmological constant divide 
w\ = —1 at a time. 


3 Fit the model parameters to Supernovae data 


The Exponential Cardassian models predict a specihc form of the Hubble 
parameter H{z) as a function of redshift z in terms of two parameters 
and n. Using the relation between cIl^z) and the comoving distance r{z) (where 
z is the redshift of light emission) 

diiz) = r{z){l +z), (18) 


and the light ray geodesic equation in a flat universe c dt = a{z) dr{z) where 
a{z) is the scale factor. 

In general, the approach towards determining the expansion history H{z) is 
to assume an arbitrary ansatz for H[z) which is not necessarily physically mo¬ 
tivated but is specially designed to give a good £t to the data for di^z). Given 
a particular cosmological model for H{z;ai,an) where Oi,..., a„ are model 
parameters, the maximum likelihood technique can be used to determine the 
best £t values of parameters as well as the goodness of the £t of the model 
to the data. The technique can be summarized as follows: The observational 
data consist of N apparent magnitudes mi{zi) and redshifts Zi with their cor¬ 
responding errors arm These errors are assumed to be gaussian and 
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uncorrelated. Each apparent magnitude m* is related to the corresponding 
luminosity distance di by 


m{z) = M + 5 logio 


dLjz) 

Mpc 


+ 25, 


(19) 


where M is the absolute magnitude. For the distant SNela, one can directly 
observe their apparent magnitude m and redshift z, because the absolute mag¬ 
nitude M of them is assumed to be constant,he. , the supernovae are standard 
candles. Obviously, the luminosity distance diiz) is the ‘meeting point’ be¬ 
tween the observed apparent magnitude m{z) and the theoretical prediction 
H{z). Usually, one dehne distance modulus fi{z) = m{z) — M and express it 
in terms of the dimensionless ‘Hubble-constant free’ luminosity distance Dl 
dehned hyD l{z) = HQdL{z)/c as 

fi{z) = 5 logio(E>L(^)) + + 0 , (20) 


where the zero offset /io depends on Hq (or h) as 


ho 


5 logio 


Mpc 


+ 25 


—5 log;^o h + 42.38. 


( 21 ) 


The theoretically predicted value D^^{z) in the context of a given model 
H{z; Oi,..., ttn) can be described by [16,17,18] 



(1 + J dz' 
0 


Hq 

H{z']ai, ...an)' 


( 22 ) 


Therefore, the best £t values for the parameters (0^,0, n) of model I are found 
by minimizing the quantity 


N 




jjX’^^^Zi) — 5 log 


Dfiz. 


10 \ 


'll f^m,0) ^) ho 


l2 


2=1 




(23) 


Since the nuisance parameter /xq is model-independent, its value from a specific 
good fit can be used as consistency test of the data [19] and one can choose a 
priori value of it (equivalently, the value of dimensionless Hubble parameter 
h) or marginalize over it thus obtaining 


B{Vtrnfl) 'fd) 


+ In 





X (^Tn,0) ^) ^(f^m,0) ^) 


C 


(24) 




Fig. 3. The observed 157 Gold SNela distance modulus along with the theoreti¬ 
cally predicted curves in model I (solid line), where we take a priori that current 
dimensionless Hubble parameter h = 0.65. 

where 


N 




- 5 logio Df{Zi-, Qrn,0, u) 


n2 


i=l 


a: 


^ - 5logio Df{zi] n) 

n) = }_^ - - 


i=l 


erf 


and 


N 




i.i V 


(25) 

(26) 


(27) 


In the latter approach, instead of minimizing x^{^m,o,n), one can minimize 
which is independent of /xq- 

We now apply the above described maximum likelihood method using Gold 
dataset which is one of the reliable published data set consisting of 157 SNela 
(TV = 157) [2]. 

In Fig. 3, we show a comparison of the observed 157 SNela distance moduli 
along with the theoretically predicted curves in model I (continuous line). It 
is obvious that the exponential Cardassian model provide a good £t to the 
observational data. 

In Fig.4, contours with 68.3%, 95.4% and 99.7% confidence level are plotted, 
in which we take a marginalization over the model-independent parameter /io- 
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Fig. 4. The 68.3%, 95.4% and 99.7% confidence contours of parameters and 
n using the Gold SNela dataset and marginalizing over the model-independent 
parameter fiQ. 

The best fit as showed in the hgure corresponds to Qm,o = 0.48 and n = 1.84, 
and the minimum value of = 173.51. Clearly, the allowed ranges of the 
parameters and n favor that there exists an effective phantom energy in 
the universe. 

As for model II, using the above marginalization method, we find the mini¬ 
mum value of cto, is 173.08 and the corresponding best-fit value of 

parameters are = 0.308, (Tq = 0.708 and 1/n = 0.0142. If we assume prior 
that flmfi = 0.3 as indicated by the observation about mass function of galax¬ 
ies, the minimum value of x^(o'o,n) is 173.178 corresponding to cxo = 0.694 
and n~^ = 0.0123. 


4 On the sound speed 


In order to guarantee the classical solution of the expansion is stable, the sound 
speed cl of classical perturbations of the total cosmological fluid around homo¬ 
geneous FRW solutions should always be greater than zero. If the expansion 
of the universe is adiabatic, the sound speed of total cosmological fluid can be 
represented by cj. = dpx/SpT- For model I, 


c 


2 

s 


-nln(G^|o)(l+ ^) 


nln(G^^o)(l + ^) ^"' + ^-1 


(28) 
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2.2 



Fig. 5. The 68%, 90% and 99% confidence level contours of parameters Clm,,o and 
n using the Gold SNela dataset and marginalizing over the model-independent 
parameter fiQ with a prior = 0-3. 

In order to assure that > 0 at any time, we should have 

nln(G-;o)(l + ;2)-3- + n-l>0 (29) 


and 


+^) - 1 < 0 . 


(30) 


From condition (29), n should be greater than 1; but condition (30) can not 
be satished for every value oi z G (—1, -l-oo). Therefore, model I should only 
be considered as an effective description for the cosmic expansion. However, 
we note that if we only require that > 0 for 2 ; > 0, i.e. , classical solution 
of the expansion is stable in the past, then we should let —1 < wta < n — 1 
which is clearly satishable. 

As for model II, the sound speed can be denoted as 


na 

( ^ ^ 

n+l n 1 

(*) 

n 

-|- n — 1 


U + aJ 

n 1 


n 

1 -1 


(31) 


It is easy to hnd that cj. is always positive if the parameters satisfy the condi¬ 
tion g” < n~^ < 1. 
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5 Discussion 


In above sections, we have investigated a new kind of Cardassian models of 
the universe,dubbed Exponential Cardassian universe. Contrary to the original 
power-law Cardassian model [9], the equation-of-state parameter w of effective 
dark fluid is dependent on time that can cross the cosmological constat divide 
tCA = — 1 from wx > — 1 to wx < — 1 as the observations mildly indicate. 

However, it is worth noting that these two models are available at much large 
scales where the matter density pm ~ Pc and evolve in the light of oc (1 -|- 
z)~'^. As discussed by Gondolo and Freese [14], the gradient of the effective 
Cardassian pressure Vpcard should be able to neglected compared to that of the 
ordinary pressure Vpm at the scales of galaxies and galaxy clusters where the 
matter density is much larger than cosmic average density. In a typical galaxy 
where pm ~ lOOpc ~ lOOpcard and Pm = Pm^"^ with velocity dispersion S = 
300 km/s, if we assume that Vpcard ^ Vpm, we should require the parameter 
n in both models considered above are greater than 3. It is remarkable that 
this value of n is compatible with the data of SNela observations. Therefore, 
these two models are both available on galactic scales. 


Observationally, the exponential Cardassian models are compatible with the 
data of SNela. The further work that should do is to check it out by other cos¬ 
mological and astrophysical observations such as CMB and large scale struc¬ 
tures (LSS) as many works (for example, [20,21]) for power-law Cardassian 
models, which is out of the scope of this paper. However, it is worth noting 
that Koivisto et al [20] pointed out that in the MP Cardassian model, the 
late integrated Sachs-Wolf effect is typically very large duo to the suppression 
of fluctuations in Cardassian fluid at late times, and is then not compati¬ 
ble with the observations of CMB unless the parameters is very close to the 
ACDM model. In the exponential Cardassian models, we assume that Cardas¬ 
sian fluctuations is induced by cold dark matter (i.e., case HI in Ref. [20]) and 
fluctuations in the cold dark matter and in the Cardassian fluid are related 
adiabatically 
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(d\ogpT\ \ 


6j' 

1 -|- Wt ’ 


(32) 


where 6 and 6t denote the fluctuations in cold dark matter and Cardassian 
fluid, respectively. According to the expressions of wt, Eq.(9) for model I and 
Eq. (15) for model H, it is not hard to find that the Cardassian fluctuations are 
both not heavily suppressed as those in MP Cardassian model which S r'U 
at late times [20]. Therefore, in the exponential Cardassian models the above 
problem is alleviated to some extent. 
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Theoretically, these models, in fact as well as power-law Cardassian models, 
are phenomenologically described in a flnid approach, therefore, we need fnr- 
ther works to find a simple fnndamental theory which can derive the form of 
g{p) presented above. 
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